The CNO Bi-Cycle 


A second way of processing !H into ^He is through a series of nuclear 
reactions involving the different species of Carbon, Nitrogen, and 
Oxygen. The principle reactions are as follows: 


120 LTH — BN +y (14.1.1) 
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WO LIH — TAN tHe (14.1.10 
Note that the net sum of these reactions is 
I2 d TI = oat He-E29 2p. 


In other words, four protons are converted into two positrons, two 
electron neutrinos, and one ^He nucleus; !?C serves only as a catalyst. 
Moreover, since the cycle can begin with any of the intermediate 
produces, it will occur if any of them exists. 





The differential equations for the full CNO bi-cyle are below. To sim- 
plify the notation, the number density symbol, N has been omitted. 
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where 7 is the lifetime of the given species against reactions with a 
proton (r = 1/AfNig), and Tg is the species lifetime against beta 
(or positron) decay. In the CNO cycle, all the reactions, with the 
exception of !"O(p, ^He)'4N, are believed to be non-resonant. 








There are several items of note in the previous equations. First, 
when !?N combines with a proton, the result is an excited compound 
nucleus, THO". This nucleus then decays into either normal !9O or 
into ^C + ^He. The probability of the former reaction depends on the 
initial energy of !9O, which, of course, depends on the energy of the 
reaction. Thus, the branching ratio, y, which appears in (14.1.11) 
and (14.1.17), is actually a mean value that is averaged over the 
Maxwellian velocity distribution. It is thus dependent (weakly) on 
temperature. Its value, y ~ 4 x 10-7, indicates that most CNO 
processing occurs through the first upper part of the bi-cycle. 


A second item to note is that !?N, !?O, and !"F can be lost either 
by beta decay or by reactions with a proton. In all cases, however, 
the timescales for beta decay are much shorter than that for nuclear 
fusion (typically seconds vs. years). Thus, unless the temperatures 
are extremely high (T > 10? K), the fusion terms in these equations 
can be dropped. 


[Note: in the late stages of stellar evolution, temperatures can reach 
10° K; when this happens, these additional reactions must be in- 
cluded. They are !?N(p, y)!4O, !?O(p, y)! F, and !"F(p,4)!5Ne. In 
each case, the product quickly decays via a positron emission: 140 — 
14N, 149p — 190 ang ?Ne — TRE — 180. If any t80 is created, it 
will then combine with a proton and produce a ^He nucleus via the 
reaction !5O(p, a) ?N ] 


A corollary to the above is that because the beta decay rates are 
so short, !1?N, !?O, and !"F quickly reach equilibrium (in a matter 
of minutes). Thus, the time derivatives of (14.1.12), (14.1.15), and 
(14.1.18) can be set to zero, and these variables can be eliminated 
from the network. The calculation of the CNO cycle then reduces to 
solving six coupled differential equations: 
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Again, these equations can be solved numerically via Runge-Kutta 
integration. However, under most circumstances, the problem can 
be simplified enormously. Examine the table (again, from Clayton 
1968), which gives the reaction lifetimes as a function of burning 
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with p in g-cm ? and 7, in years. For the higher temperatures 
(which correspond to higher mass stars), the timescales listed are 
all much shorter than the timescale for the depletion of hydrogen, 
or for changes in the stellar temperature and density. Thus, while 
the CNO isotopes are achieving equilibrium, it is fair to assume that 
the hydrogen abundance and main sequence structure of the star are 


3 


constant. This is important, since it allows us to fix the values of 7 
in equations (14.1.26) - (14.1.31), and linearize the equations. 





To see this, let the subscripts 1 through 6 represent !?C, !?C, 14N, 
15N, 160. and tO, and let y be the number the density of the species. 
In differential form for back-differencing, equations (14.1.20) through 
(14.1.25) are 
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These equations can now be re-written into a more useful form: 
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which, in matrix form is 
X 0 0 X 0 U V, Ui 
X X 0 0 0 0 y^ y2 
0 X X 0 0 X ys | |% 
0 0 X X 0 0 yl | ya 
0 0 0 X X 0 Ys Ys 
0 0 0 0 X X Uo UG 


where the X's represent the terms of (14.1.32) - (14.1.37). The up- 
dated CNO abundances can therefore be found via a simple 6 x 6 
matrix inversion. 





The time it takes the CNO cycle to reach equilibrium is a critical 
function of temperature. However, two properties of the chain should 
be noted. First, the CNO equations really consist of two separate 
cycles: a CN cycle and an ON cycle. If you examine the lifetimes of 
the CN reactions, you can see that this cycle proceeds fairly rapidly, 
with the timescale for ZN fusion (r ~ 10° years at Tg ~ 20 K) 
being the limiting factor. Since the branching ratio out of this cycle 
is very small (y ~ 4 x 1074), ~ 10? cycles are required before a 
significant number of CN nuclei are lost to the ON chain. Meanwhile, 
the timescale for ON nuclei to filter into the CN cycle is similarly long, 
T ~ 108 years (the !6O reaction timescale for Tg ~ 20 K). As a result, 
the CN cycle falls into quasi-equilibrium fairly quickly, and is largely 
decoupled from the slower ON cycle. 











The second thing to note is that in high mass stars (i.e., stars with 
high central temperatures), the reaction timescales are short, and 
the cycle achieves equilibrium fairly quickly. Thus, for purposes of 
estimating energy generation, the equilibrium abundances of the el- 
ements can (most times) be used. On the other hand, in low mass 
stars, where CNO equilibrium may never be achieved, the proton- 
proton chain dominates energy production. Thus, small errors in the 
CNO reaction are acceptable. 








Of course, the operation of the CNO cycle means that the differential 
equations for 'H and ^He must be modified to include the additional 
destruction and creation mechanisms. 


Table 5-3 Dependence of log (rp 4/100) on temperaturet 


Temper- Reaction] 
00, 
T, | CU(p,y)N' Cip, y) N Nip, yO! N'(pa)Ci: 10y = O**(p,y}F 7 O27 (pce) NM 
5 16.32 15.73 19.79 15.53 4.649 22.95 21.92 
6 14.32 13.73 17.57 13.29 4.598 20.51 20.02 
7 12.72 12.13 15.79 11,50 4.551 18.56 . 18, 26 
8 11.41 10.81 14.32 10.03 4,508 16.95 16,50 
9 10.29 9.69 13.08 8.78 4,468 15.59 15.10 
10 9.33 8.73 12.02 7.70 4.431 14.42 14.05 
11 8.50 7.90 11.09 6.76 4.396 13.39 13.15 
12 7.75 7.15 10.26 5.93 4.363 12.49 12.38 
13 7.09 6.49 9.52 5.18 4.332 11.68 11.68 
14 6.49 5.89 8,86 4.51 4.303 10.95 11.02 
15 5.95 5.35 8.26 3.90 4.275 10.29 10.32 
16 5.45 4.85 7.71 3.34 4,248 9.68 9.55 
17 5.00 4.39 7.20 2.83 4.223 9.13 8.70 
18 4.58 3.97 6.73 2.35 4.198 8.61 7.86 
19 4.18 3,58 6.30 1.91 4.175 8.14 7.01 
20 3.82 3.21 5.89 1.50 4.152 7.69 6.18 
22 3.16 2.95 5.16 0.7 4.110 6.89 4.78 
24 2.57 1.97 4.51 0.09 4.071 6.18 3.63 
25 2.30 1,70 4.21 —0.21 4.052 5.85 3.10 
26 2.05 1.44 3.93 —0.50 4.034 5.54 2.62 
28 1.58 0.97 3.41 —1.03 4.000 4.97 1.75 
30 1.15 0.54 2.93 —1.31 3.967 4.45 1.05 
35 0.23 —0.38 1.91 —2.55 3.803 3.33 —0.42 
40 —0.53 —1.14 1.07 — 83.42 3.829 2.41 —1.50 
45 —1.18 —1.78 0.36 —4.14 3.771 1.64 —2.33 
50 —1.73 — 2.43 —0.25 —4.77 3.719 0.97 —2.99 
55 —2.21 —2.82 —0.78 —5.32 3.673 0.39 —3.53 
60 —2.64 —3.24 —1.25 —5.81 3.630 —0.12 —3.97 
65 —3.02 —3.63 —1,67 —6.24 3.590 —0.58 —4.33 
70 —3.37 —3.97 —2.05 —6.63 3.554 —~ 0.99 —4.65 
75 — 3.68 —4.28 —2.39 —6.99 3.521 —1.37 —4.91 
80 —3.97 —4.57 —2.71 —7.32 3.489 —1.71 —5.14 
85 —4.23 —4.83 —2.99 —1.62 3.460 —2,02 —5.35 
90 — 4.48 —5.08 ~3.26 —7.90 3.433 —2.31 —5.52 
95 —4.70 —5.30 —3.51 -8.15 3.407 —2.98 —5.68 
100 —4.91 —5.51 —3.74 —8.39 3.383 —2.83 —5.82 


ee EN RS 
7 Adapted from G. R. Caughlan and W. A. Fowler, Astrophys, J., 136:453 (1962). By permis- 
sion of The University of Chicago Press, Copyright 1962 by The University of Chicago. 


į The lifetimes against protons are expressed in years, and the density p is in grams per cubic 
centimeter. 


Energy Generation in the CNO Cycle 


The energy produced by the CNO cycle is simply the sum from the 
individual reactions in the chain. However, we can approximate the 
energy generation in equilibrium by neglecting the ON chain, and 
considering that the *cycle time" for CN is dominated by the slowest 
reactions. Thus, we can approximate the chain with three steps: 


*C(p, y) ?5N(8*, v)P?C 


I3O(p, y)14N 
AN (p, 4)90(8* , v) *N(p,! He)?C 


The first reaction produces 1.944 + 2.221 — 0.710 = 3.455 MeV per 
reaction; the second generates 7.550 MeV per reaction, and the third 
14.019 MeV, per reaction. The total energy for the chain in equilib- 
rium is therefore 


N(**C NTL NGAN 
T12 T13 T14 
x 1079 ergs s! cm? (14.2.1) 


or, since the abundance of the individual isotope is simply propor- 
tional to its lifetime. 
N(CN 
pec © 40.08 x 109 R M ergs s ! cm? (14.2.2) 
Ti»er Tae r5 

Now recall that these are non-resonant reactions. The nuclear life- 
times will depend on the reaction rates, A, which has a temperature 
dependence of 
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Note that this dependence is proportional to Z?/?. 'Thus, we can 
expect that the CNO cycle's energy generation will be ^ 5 times 
more sensitive to temperature than the proton-proton. 


A reasonable approximation to the CNO cycle’s rate of energy pro- 
duction in equilibrium is 


eng = 8.67 x 10?” gia Xono X p Tg O e 15228/75^* ergs st em 3 
with 
9141 = 1 + 0.0027 T$? — 0.00778 T2/? — 0.000149 T; 


As Tg varies from 10 to 50, v goes from 23 to 13. This compares to 
V ~ 4 for the proton-proton chain. 


N24 
N 
N x 104 
0.1 C2? x 10 

o 

LS) 

d 

a 

v 

c 

E ee 

e" 

7 

E ) 

3 

uc 

tJ 

0.01 
e104 
Ò. 16 
e 
O27 x 10 
0.001 
0 20 40 60 80 100 


Temperature, 105*K 


CNO Equilibrium Abundances 


Ratio of equilibrium abundances 





60 
Temperature, 10°°K 


CNO Equilibrium Isotope Ratios 


Relative abundance 


0.001 


€ 
pXyu* 


log ( 


; (PX 
Time, years T(P Ando 

100 
b N m m * ooa n o LJ ~ L3 L.] [LJ 
o o o o oooo0 o o o o o o 
Laa — — — Lu m ~ — m — E rd 
x x x x X Xx x x x x x x x 
Lor e — m vd Oe H -= m ~ [12] — r 





9 
3 


o 
e 
— 


0.01 


Protons consumed per initial nucleus 


CNO approach to equilibrium starting with 50% !^C, 50% !9O. 


8 


CNO bi-cycle 


PP chain 





10 15 20 25 


Ts 


30 35 40 45 


CNO vs. PP energy generation for Xow/ Xg = 0.02. 


The Neon-Sodium Cycle 


The third way of turning !H into ^He is through the Ne-Na cycle. 
This cycle only occurs in extremely massive main-sequence stars (at 
temperatures of Tg ~ 90 K) and cannot compete energetically with 
the CNO cycle. However, it is important for changing the ratios 
of these isotopes, especially 77Ne which can be an important source 
of neutrons for heavy element nucleosynthesis. The reactions of the 
Ne-Na cycle are 

















Ne + 1H — Na +y (14.3.1) 

21Na — ?!Ne + e* + ve (14.3.2) 
^INe + 1H — Na +y (14.3.3) 

22Na — ?2Ne + et + v, (14.3.4) 
22Ne + 1H — Na +y (14.3.5) 
23Na + +H — ??Ne + *He (14.3.6) 

or 

23Na + 1H — ?*Mg +y (14.3.7) 


Compared to the strong-force interactions, the beta-decay reactions 
occur almost instanteously, so (14.3.2) and (14.3.4) can immediately 
be merged into (14.3.1) and (14.3.3). Also, note that the Ne-Na cycle 
is not really a complete cycle, since some ??Na nuclei will be fused to 
^^Mg. The split between these two branches is about 50-50. Thus, 
after a few cycles, Ne and Na will become depleted, in favor of the 
more-inert isotope ?^Mg. 





